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A new generation of local three-body potentials providing an excellent description of the properties 
of light nuclei, as well as of the neutron-deuteron doublet scattering length, has been recently derived. 
We have performed a comparative analysis of the equations of state of both pure neutron matter 
and symmetric nuclear matter obtained using these models of three-nucleon forces. None of the 
considered potentials simultaneously explains the empirical equilibrium density and binding energy 
of symmetric nuclear matter. However, two of them provide reasonable values of the saturation 
density. The ambiguity concerning the treatment of the contact term of the chiral inspired potentials 
is discussed. 

PACS numbers: 21.30.Fe, 21.45.Ff, 21.65.-f 
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I. INTRODUCTION 

The definition of the potential describing three-nucleon 
interactions is a central issue of nuclear many-body the- 
ory. Three-nucleon forces (TNF) are long known to pro- 
vide a sizable contribution to the energies of the ground 
and low-lying excited states of light-nuclei, and play a 
critical role in determining the equilibrium properties of 
symmetric nuclear matter. In addition, their effect is ex- 
pected to become large, or even dominant, in high density 
neutron matter, the understanding of which is required 
for the theoretical description of compact stars. 

Ab initio nuclear many-body approaches are based on 
the premise that the dynamics can be modeled studying 
exactly solvable systems, having mass number A < 5. 
This is a most important feature since, due to the com- 
plexity of strong interactions and to the prohibitive diffi- 
culties associated with the solution of the quantum me- 
chanical many-body problem, theoretical calculations of 
nuclear observables generally involve a number of approx- 
imations. Hence, models of nuclear dynamics extracted 
from analyses of the properties of complex nuclei are 
plagued by the systematic uncertainty associated with 
the use of a specific approximation scheme. 

Highly realistic two-nucleon potentials, either purely 
phenomenological |lj-y| or based on chiral perturbation 
theory (ChPT) 0, m|, have been obtained from accurate 
fits of the properties of the bound and scattering states of 
the two-nucleon system [vl-fisj. Unfortunately, however, 
the extension to the case of the three-nucleon potential 
is not straightforward. Phenomenological models, such 
as the Urbana IX (UIX) potential, that reproduce the 
observed binding energy of ^H by construction, fail to 
explain the measured nd doublet scattering length, '^Und 
[14| . as well as the proton analyzing power in p-^Hc scat- 
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In recent years, the scheme based on ChPT has been 
extensively emp loyed to obtain three-nucleon potential 
models [la, ll3 ■ The main advantage of this approach 
is the possibility of treating the nucleon-nucleon (NN) 
potential and the TNF in a more consistent fashion, as 
the parameters ci, C3 and C4, fixed by NN and ttN data, 
are also used in the definition of the TNF. In fact, the 
next-to-next-to-leading-order (NNLO) three-nucleon in- 
teraction only involves two parameters, namely cd and 
ce, that do not appear in the NN potential and have to 
be determined fitting low-energy three-nucleon (NNN) 
observables. Unfortunately, however, ttA^ and NN data 
still leave some uncertainties on the Ci's. that can not be 
completely determined by NNN observables. 

A comprehensive comparison between purely phe- 
nomenological and chiral inspired TNF, which must nec- 
essarily involve the analysis of both pure neutron matter 
(PNM) and symmetric nuclear matter (SNM), is made 
difficult by the fact that chiral TNF are derived in mo- 
mentum space, while many theoretical formalisms are 
based on the coordinate space representation. 

The local, coordinate space, form of the chiral NNLO 
three nucleon potential, hereafter referred to as NNLOL, 
can be found in Ref. [I3. However, establishing a con- 
nection between momentum and coordinate space repre- 
sentations involves some subtleties. 

The authors of Ref. ^ have shown that the NNLO 
(momentum space) three body potential obtained from 
the chiral Lagrangian, when operating on a antisymmet- 
ric wave function, gives rise to contributions that are not 
all independent of one another. To obtain a local po- 
tential in coordinate space one has to regularize using 
the momenta transferred among the nucleons. This reg- 
ularization procedure makes all the terms of the chiral 
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potential independent, so that, in principle, all of them 
have to be taken into account. The potential would oth- 
erwise be somewhat inconsistent, as it becomes apparent 
in nuclear matter calculations, which involve larger mo- 
menta. 

A comparative study of different three- nucleon local in- 
teractions (Urbana UIX (UIX), chiral inspired revision of 
Tucson-Melbourne (TM') and chiral NNLOL three body 
potential), used in conjunction with the local Argonne 
Vis NN potential, has been recently performed [l^. The 
authors of Ref. [l9| used the hyperspherical harmonics 
formalism to compute the binding energies of ■^H and 
"^He, as well as the nd doublet scattering length, and 
found that the three body potentials do not simultane- 
ously reproduce these quantities. Selecting different sets 
of parameters for each TNF they were able to obtain 
results compatible with experimental data, although a 
unique parametrization for each potential has not been 
found. This problem is a consequence of the fact that the 
three low-energy observables considered are not enough 
to completely fix the set of parameters entering the defi- 
nition of the potentials. 

The work described in this paper is aimed at testing 
the different parametrization of the potentials in nuclear 
matter. In the case of SNM, a realistic Equation of State 
(EoS) is constrained by the available empirical informa- 
tion on saturation density, po, binding energy per nucleon 
at equilibrium, i?o, and compressibility, K. Furthermore, 
the recent observation of a neutron star of about two solar 
masses [2^ puts a constraint on the stiffness of the EoS 
of beta-stable matter, closely related to that of PNM. 

Nuclear matter calculations are carried out using a va- 
riety of many-body approaches. The scheme referred to 
as FHNC/SOC; based on correlated basis functions and 
the cluster expansion technique, has been first used to 
perform accurate nuclear matter calculations with realis- 
tic three body potentials in Ref. [2l|. This analysis in- 
cluded early versions of both the Urbana (UIV, UV) and 
Tucson Melbourne (TM) three body interactions with 
the set of parameters reported in Ref. [1^. The results 
indicate that the UV model, the only one featuring a phe- 
nomenological repulsive term, provides a reasonable nu- 
clear matter saturation density, while the UIV and TM 
potentials fail to predict saturation. In addition, none 
of the considered models yields reasonable values of the 
SNM binding energy and compressibility. 

The findings of Ref. [2lj are similar to those obtained in 
Ref. [l^ , whose authors took into account additional dia- 
grams of the cluster expansion and used the UVII model. 
The state-of-the-art variational calculations discussed in 
Ref. [i^l, carried out using the Argonne wis and UIX 
[25} potentials, also sizably underbinds SNM. While the 
authors of Ref. 2J] ascribed this discrepancy to defi- 
ciencies of the variational wave function, the analysis of 
Refs. [2^ [13 suggest that this problem can be largely 
due to the uncertainties associated with the description 
of three-nucleon interactions, whose contribution turns 
out to be significant. 



Momentum space chiral three-body interaction have 
been also employed in nuclear matter [28l - [30| . In these 
studies, the NNNLO chiral two-body potential has been 
evolved to low momentum interaction Viowk, suitable 
for standard perturbation theory in the Fermi gas basis. 
The results, showing that the TNF is essential to obtain 
saturation and realistic equilibrium properties of SNM 
[28. 29], exhibit a sizable cutoff dependence. At densi- 
ties around the saturation point this effect is ~ 4MeV. 
In addition, different values of the constants Ci lead to 
different Equations of State for SNM [H and PNM 

The main features of the chiral inspired TNF are briefly 
reviewed in Section|TIl while in Section HIIl we analyze the 
coordinate space form of the TNF derived in Ref. [l^ 
and discuss several issues related to the calculation of 
their contributions in nuclear matter, both within the 
FHNC/SOC and the Auxihary Field Diffusion Monte 
Carlo (AFDMC) approaches. The numerical results, in- 
cluding the EoS of PNM and SNM are discussed in Sec- 
tion |lVl Finally, in Section |V] we summarize our findings 
and state the conclusions. 



II. CHIRAL INSPIRED MODELS OF THREE 
NUCLEON FORCES 

In a chiral theory without A degrees of freedom, the 
first nonvanishing three-nuclon interactions appear at 
NNLO in the Weinberg power counting scheme [3l|, [s^ . 
The interaction is described by three different physical 
mechanisms, corresponding to three different topologies 
of Feynman diagrams, drawn in Fig. [T] [l6j . The first 
two diagrams correspond to two-pion exchange (TPE) 
and one-pion exchange (OPE) with the pion emitted (or 
absorbed) by a contact NN interaction. The third dia- 
gram represents a contact three-nucleon interaction. 

The full expression for the TNF is obtained by sum- 
ming all possible permutations of the three nucleons. For 
this kind of potential, it turns out that there are only 
three independent cyclic permutations, i.e. 

V^{1, 2, 3) = y>^(l : 2, 3) V^^(2 : 1, 3) 

+ ^^^(3:1,2) (1) 

The Feynman diagrams of Fig. [T]refer to the permutation 
(3 : 12), that can be written as 

yx(3 : 12) = 01^1(3 : 12) + €3^^3(3 : 12) + 04^4(3 : 12) 
+ cdVd {3 : 12) + ceVe{3 : 12) . (2) 

The first three terms Vi , V3 and V4 come from the TPE 
diagram and are related to ttN scattering. In particular, 
Vi describes the S'-wave contribution, while V3 and V4 are 
associated with the P-wave. The other terms, Vjj and 
Ve, are the OPE and contact contributions, respectively. 
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Figure 1. TPE, OPE and NNN contact interactions of the 
chiral three-body force at NNLO. 



Their momentum space expressions are [16| 
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where Ti and are the Pauli matrices describing the 
spin and the isospin of particle i. With (Ty- and r^j we 
denote the scalar product cr j • crj and • Tj , respectively. 
The strengths of the TPE, OPE and contact terms, Vq, 
Vff and Vff are given by 
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= 92.4 MeV is the weak pion decay constant and 
is the chiral symmetry-breaking scale, of the order of the 
p meson mass. 

The low energy constants (LEG) ci, C2 and C3 also ap- 
pear in the sub-leading two-pion exchan ge term of the 
chiral NN potential and are fixed by TriV [331. [3^ and/or 
NN Q data. The parameters cjj and ce are specific 
to the three-nuclcon interaction and have to be fixed us- 
ing NNN low energy observables, such as the '^H binding 
energy and the nd doublet scattering length ^a„d p^ . 

The many-body methods employed in our work, 
namely FHNC/SOC and AFDMC, require a local ex- 
pression of the three-body potential in coordinate space, 
that can be obtained performing the Fourier transform 

m 
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where the cutoff functions Ff^, defined as 
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can depend on the momenta transferred among the nu- 
cleons, g,, only. This feature has important consequences 
for the OPE and contact terms, that will be discussed at 
a later stage. 

The cutoff A in the previous equation, while not being 
required to be the same as A^, is of the same order of 
magnitude. Choosing the fourth power of the momen- 
tum in Eq. ^ is therefore convenient, as the regulator 
generates powers of q/ K which arc beyond NNLO in the 
chiral expansion. 

The Fourier transform can be readily computed, and 
provides the following coordinate-space representation of 
the chiral three-body potential: 

Vi{2, : 12) = Wq ti2((Ti • 7%3)(cr2 ■ ^23)2/(^13)2/(^23) 

V3(3 : 12) = VKo Ti2[(Ti22/(ri3)2/(?'23) 

+ (CTI • ?^23)('7'2 • r2z)t{r23.)y{ri3.) 

+ (cti • f^i3)(o-2 • riz)t{riz)y{r2z) 

+ (na • ■r2z){<Ji ■ ri3)(o-2 • ■r2z)t{riz)t{r2z)] 

14(3 : 12) = Wn (rg • ti x X2)[(cr3 • 0-2 x (Ti)y{ri'i)y{r2z) 
+ {cr-A ■ r23 X cri)((T2 • r23)t(r23)2/(?'i3) 
+ (0-2 ■ fi3 X cr3)((Ti • ri3)t(ri3)2/(r23) 
+ (era ■ r23 x fi3)((Ti • fi3)(cr2 • r23)t{ri3)t{r23)] 

12) = Wo^Ti2[(Ti22/(r23)zo('"l3) 

+ (ci • ■r23)i(T2 ■ r23)tir23)zo{ri3) 
+ cri2y{ri3)M'i^23) 
+ {0-2 ■ ri3){cTi ■ ri3)t(ri3)zo(r23)] 

12) = Win2ZQ{ri3)zQ{r23) , 



(7) 



where Wq, VFjf and are obtained multiplying the 
where = 1-29 is the axial-vector coupling constant, corresponding Vq, Vff and Vjj^ by a factor m^/(47r)^. 
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The radial functions appearing in the above equations 
are defined as 



y{r) 

t{r) 



z[{r) 



r 



r 



-y'{r) 



(8) 



while Zm proportional to Z„ introduced in Ref. [22 
given by 
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with jo{x) = sin(a;)/a;. Note that, due to the form of the 
cutoff function of Eq. the radial functions are not 
known in analytic form, and must be obtained from a 
numerical integration. 

Recently, the authors of Ref. [3 have studied the low 
energy NNN observables using the hyperspherical har- 
monics formalism and a nuclear hamiltonian including 
the NNLOL potential and the Argonne vis Q two-body 
interaction. This mixed approach requires a fit of all 
the LEG appearing in the chiral three-body interaction, 
not CD and ce only. The best fit parameters for the 
■^H and "^He binding energies and for the nd scattering 
length, '^Und, are listed in Table ID For all the differ- 
ent parametrizations, denoted by NNLOL^, ci and A^^, 
have been fixed to their original values 0.00081 Me 
and 700 MeV, respectively [l6|. The momentum cutoff 
of Eq. (O has been set to 500 MeV. 

Table I. Parameters of the NNLOL interactions of Ref. p^ . 



Potential C3 (MeV"^) 


C4 (MeV"') 


CD 


Ce 


NNLOLi -0.00448 


-0.001963 


-0.5 


0.100 


NNLOL2 -0.00448 


-0.002044 


-1.0 


0.000 


NNLOL3 -0.00480 


-0.002017 


-1.0 


-0.030 


NNLOL4 -0.00544 


-0.004860 


-2.0 


-0.500 



As noticed in Ref. [37[ , despite the different underlying 
physical mechanisms, both TM and UIX three-nucleon 
interactions can be written as a sum of terms of the 
same form as those appearing in Eq. d?]). The differ- 
ences among NNLOL, TM and UIX lie in the constants 
and in the radial functions. 

The TM' potential only involves the Vi, V3 and V4 
contributions [35|. The cutoff function for this potential 
is not the same as in Eq. ((H]), but 
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The above form allows for the analytical integration of 



Eq. ([9]), yielding the radial functions 
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r(A— ttItt ) 



-rA 



- 6 + 2(3 + 'im^r + mlr^)e''^'^^"'-'^ 



+ m^r2(l + rA) -rA[6 + rA(3 + rA)] . (11) 

The TM' potential corresponds to the following choice 
of the strength constants (compare to Eq. (O) 
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a, b and c being the parameters entering the definition of 
the TM' potential [3|. The authors of Ref.[ii| have de- 
termined the parameters of the TM' potential fitting the 
same set of low energy NNN observables employed for the 
NNLOL potential. In order to get a better description of 
the experimental data, they introduced a repulsive three- 
nucleon contact term, similar to the chiral Ve, but with 
T12 omitted 



Ve{3: 12) = W^zo{r 13) zo{r23), 



where 



2 9m^ 



VSTrmAT/ A^ 



(14) 



(15) 



The corresponding radial function can be computed an- 
alytically from Eq. ([9]) 



,-rA 



Mr) = 



SttA 



(ml - Ar 



(16) 



As in the original paper [23l, in Ref. [13] the value 
of the pion-nucleon coupling constant is set to = 
179.7 MeV, the pion mass is 111^^ ~ 139.6 MeV and the nu- 
cleon mass is defined through the ratio mfq jm-^ — 6.726. 
The symmetry breaking scale A^ of Eq. has the same 
value, 700 MeV, used for the NNLOL potential. 

The parameters of the TM' potentials, TM,^, that ac- 
cording to Ref.[l^ reproduce the binding energies of "^H 
and ^He and ^a„d, are listed in Table [III It turns out 
that Fl, gives a very small contribution to the low en- 
ergy NNN observables. Therefore, the parameter a has 
been kept to its original value — 0.87to~^. 

The Fujita Miyazawa term [s^ of the UIX potential 
[2^, V'^'^ ^ describing the process whereby two pions are 
exchanged among nucleons and a A resonance is excited 
in the intermediate state, is conveniently written as 



1/2- (3 : 12) = A2^{Ai3,A23}{ri3,r23} 

-I- C27r [Ai3 , X23] [ri3 , r23] , 



(17) 
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Table II. Parameters of the TM' potential reproducing low 
energy the NNN experimental data with a — — 0.87m^^ [l9l |. 



Potential b{m^^) 


d(m-3) 


ce 


A(m^) 


TM'i -8.256 


-4.690 


1.0 


4.0 


TMa -3.870 


-3.375 


1.6 


4.8 


TM3 -2.064 


-2.279 


2.0 


5.6 



was obtained from the fit of the observed binding energy 
of ■^H, the strength Uq, was indeed adjusted to reproduce 
the empirical saturation density of SNM, po = 0.16 fm~'^. 

In Ref. [l^ it has been found that the original 
parametrization of the UIX potential underestimates 
^Und and slightly overbinds of ^He. 



III. THREE NUCLEON POTENTIALS IN 
NUCLEAR MATTER 



where 



and 



Sij = 3(fy • cri){fi 



(18) 



(19) 



is the tensor operator. The radial functions associated 
with the spin and tensor components read 



The investigation of uniform nuclear matter may shed 
light on both the nature and the parametrization of the 
TNF, although the quantitative description of this sys- 
tem can not be achieved within a mere generalization of 
the approaches developed for light nuclei. In this Sec- 
tion, we analyze the structure of the contact term of the 
NNLOL potential of Ref. [l^ and discuss the calculation 
of the TNF contribution to nuclear matter energy. 



Y(x) = ^Yix) 

X 



(20) 



T{x)=(i + - + -)y{x)^t{x) (21) 
and the ^(x) are short-range cutoff functions defined as 

iYix)=^T{x) = l-e-'=^\ (22) 

In the original derivation of the UIX potential the ratio 
C27r/^27r was fixed to 1/4 and the cutoff parameter was 
c = 2.1 fm^^, the same value as in the cutoff functions of 
the one-pion exchange term of the Argonne wis two-body 
potential. 

It can be shown that the anticommutator and commu- 
tator terms correspond to V3 and V4 of Eq. ([7|) , provided 
the following relations between the constants 



bWo = 4A2. 
dWo = 

and the radial functions 

f Y{r) =y{r) + ^-^t{r) 
\ T{r) = rlt{r) 

are satisfied. 

The repulsive term of the UIX potential 

1/^(3 : 12) = UoT\m^n3)T\m^r23) 



(23) 



(24) 



(25) 



is equivalent to the Ve term appearing in the TM' poten- 
tial and (aside from the T12 factor) in the NNLOL chiral 
potential if the following relations hold 



T^im^r) = zo(r) , Uo = ceWo^ 



(26) 



The UIX potential was not designed to reproduce low 
energy NNN obscrvables only. While the parameter A2tt 



NNLOL contact term issue 



While the NNLOL chiral interactions provide a fully 
consistent description of the binding energies of '^H and 
"'He, as well as of the scattering length ^a„d, some ambi- 
guities emerge when these interactions are used to calcu- 
late the nuclear matter EoS. 

For our purposes, it is convenient to rewrite the 
NNLOL chiral contact term of Eq, ([7]) in the form 



yi(3 : 12) ^ Vo''n2Zo{ri3)Z„{r23) . 



(27) 



where the superscript r has a meaning that will be soon 
clarified. The radial function Zo{r) = m^/(47r)2:o(r) ap- 
proaches the Dirac 5-function in the limit of infinite cut- 
off. Strictly speaking, the local version of Ve is a genuine 
"contact term" in this limit only, while for finite values 
of the cutoff it acquires a finite range. 

In addition to of Eq. ((27|) . the chiral expansion 
leads to the appearance of five spin-isopin structures in 
the contact term. For example, the scalar contribution is 



l/i(3:12) = yo^Zo(ri3)Zo(r23). 



(28) 



Within this context, the superscripts r and / identify the 
T12 and scalar contact terms. 

In Ref. [1^ it has been shown that, once the sum 
over all cyclic permutation is performed, all contribu- 
tions to the product between the potential and the anti- 
symmetrization operator ^123 have the same spin-isospin 
structure. Therefore it is convenient to take into account 
just one of the contact terms. This result was obtained 
in momentum space, without the cutoff functions F\. As 
a consequence, in coordinate space it only holds true in 
the limit of infinite cutoff. In particular, for Vj^{3 : 12) 
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and V^(3 : 12), it turns out that 



implying that in PNM 



cycl cycl 

(29) 

making this two terms equivalent. The limit of infinite 
cutoff is crucial, because the radial part of the exchange 
operator, when multiplied by the Dirac (5-functions, is 
nothing but the identity 



1 



(30) 



After the regularization, i.e. with the (5-function replaced 
by Zq, the proof is spoiled and the six different structures 
are no longer equivalent. 

In PNM contact terms involving three or more neu- 
trons vanish because of Pauli principle. On the other 
hand, the expectation value of the contact terms of the 
NNLOL potential can be different from zero. 

Let us assume that reproducing the binding energies 
of light nuclei and ^Und require a repulsive Ve ■ Then, 
one has to choose either c^^ < or cjj > 0. In PNM, as 



(n 



2] PNM 



= 1, 



(31) 



it turns out that is attractive and repulsive. This 
means that fitting the binding energies and the n ^ d 
scattering length with either or alone leads to an 
ambiguity in the expectation value of the potential. 
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Figure 2. Radial dependence of the Function Zo{r), appearing 
in EQ. (|27p . plotted for different values of the cutoff A. 



By expanding the cutoff function 



(1/i^"(3 : 12))pjvA/ 



(34) 



From the above equation it becomes apparent that the 
expectation value of the three-nucleon potential, as well 
as its sign ambiguity, is nothing but a a cutoff effect. 
Hemce, it should be regarded as a theoretical uncertainty. 
Note that, since A^^, ~ A, then {Ve)pnm is of the same 
order of the next term in chiral expansion. 

To clarify this issue, let us consider a simple system: 
a Fermi gas of neutrons, in which correlations among 
particles are not present. The expectation value of the 
contact interaction reads 



I PNM 



= y^o^/ c^^7-i2dVi3Zo(ri2)Zo(ri3)x 

^ t{ri3? i{r23? 

2 2 2 

2 J' 



(35) 



where A is the number of neutrons. The factor 1/2 in- 
cludes the 1/3! arising from the unrestricted sum over 
particle indices 123, multiplied by a factor 3 from the 
cyclic permutations of the potential, all giving the same 
contribution. The Slater function i{rij), for a system of 
fcrmions with degeneracy d is given by 



|k|<fcj=. 

sm{kFrij) — kprij coslkprij) 



(36) 



where kp = {^^^Y^'^ is the Fermi momentum. It can be 
easily seen that, if V^'^(l : 23) cx 5(?'i2)^(?'i3), then 



\^ El PNM 
A 



= 0. 



(37) 



Consider now a Fermi gas with equal numbers of pro- 
tons and neutrons, where 



SNM 
A 



^y^f y d^n2d\i3Zo{ri2)Zo{n3)x 
-liir23Y + liin2)i{n3)iir23)) (38) 



A* VA8 



one finds 

V^{3 : 12) = y„^Ti2 [s{n3)S{r23) + o[^) 
Vii3:12) = V,^[5ir^3)6ir23) + 0[^) 



(32) 



(33) 



and 



\^E/SNM 

A 



Y^o^ I d^n2d^n3Zoiri2)Zo(,n3)x 



eiri2Y eiruY l{r23? 



1 - 

4 4 

l{TX2)i{rX3)i{r23)\ 



(39) 
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Table III. Cutoff dependence of tfie expectation values of the 
three body contact term of the NNLOL potential in nonin- 
teracting PNM. 



A(MeV) (Vi-"i^)gg^,M(MeV) 



300 


9.15 


400 


5.95 


500 


3.60 


600 


2.15 


700 


1.30 


800 


0.81 


oo 






Table IV. Same as in Table |llll but for SNM. 



A(MeV) (V-^")g^MM(MeV) (MeV) 



300 


-2.61 


10.21 


400 


-3.61 


8.15 


500 


-4.37 


6.93 


600 


-4.87 


6.30 


700 


-5.15 


5.98 


800 


-5.30 


5.81 


oo 


-5.55 


5.55 



In the limit of infinite cutoff the above equations imply 
A " "16^ ^° 

As expected from Eq. ([29]) . the two contributions have 
opposite sign. 

We have computed the expectation values of Eqs. ([55]) . 
(pSj) and (f39|) for different values of the cutoff A and den- 
sity p — 0.16 fm~^. The results listed in Table Hill show 
that for PNM the larger the cutoff the smaller the expec- 
tation value of the three nucleon contact term. Note that 
for A = 500 MeV, the expectation value is still sizably 
different from the asymptotic limit. 

As far as SNM is concerned (see Table HV]). as the cutoff 
increases the possible choices of the three nucleon contact 
term tend to the asymptotic values of Eq. (|40|). As in the 
case of PNM, the results corresponding to A = 500 MeV, 
are significantly different from the asymptotic values. 

We emphasize that the parameter ce has not been in- 
cluded in this analysis, even though it is itself cutoff de- 
pendent. Unfortunately, the authors of Ref. [l^ kept A 
fixed to 500 MeV. Had this not been the case, their fit 
to the experimental data would have resulted in a set of 
different constants ce, corresponding to different values 
of A. It would have been interesting to extrapolate the 
expectation value of Ve to the limit of infinite A, where 



the cutoff effects associated with the regularization pro- 
cedure are expected to vanish. 



B. FHNC/SOC calculations 

The diagrams involved in the FHNC/SOC calculation 
of the expectation values of the V"^^ and terms of 
the UIX potential are depicted in Figs. [3] and |4l re- 
spectively. The thick lines represent the potential, while 
dashed and wavy lines correspond to generalized scalar 
and opcratorial correlations, denoted by Z'^ and in 
Ref. \M. Double wavy lines represent Single Operator 
Rings (SOR), while vertex corrections, although included 
in the calculations, are not shown. The definitions of all 
these quantities can be found in Refs.fsollioj. 

Note that, because of the symmetry properties of the 
wave function, we can restrict our analysis to the per- 
mutation (3 : 12). Taking into account the other per- 
mutations results in the appearance of a multiplicative 
factor. 

The computation of all of diagrams (3. a), (3.b) and 
(3.c) and all diagrams of Fig. |3]is outlined in Ref. [2l|, 
while the contribution of digram (3.d), involving three 
non central correlations was first taken into account by 
the authors of Ref. [2^ . 



3 3 




Figure 3. Cluster diagrams contributing to the expectation 
value of V^"". 

Using the relations for the constants and the radial 
functions given in Eqs. ((23)) and ([24]), the computation 
of the diagrams of Fig. [3] with the V3 and V4 terms of 
both the TM' and NNLOL potentials is the same as that 
of 1/2^ reported in Ref. [2l|. 
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where 






Figure 4. Same as in Fig. [3l but for V^. 



Thanks to the identity 

(ci • ri3)(cr2 • f23)(ri3 • ^23) = 
Y^{ci3 + Sis,a23 + S23} , 



(41) 



the term Vi of Eq. ([7]), appearing in both the TM' and 
the NNLOL potentials, can be written in the form 



^1(3 : 12) 



^ Wo 

36 
ri3r23 



{■''13,'''23}{0'13 + <S'i3,(T23 + 'S'23}x 

y{ri3)y{r23) ■ (42) 



{fl3 ■ r23) 

Aside from the radial function, Vi is completely equiva- 
lent to V3, the anticommutator term of the UIX poten- 
tial. Therefore, we were allowed to use again the results 
of Ref. lU. 

Furthermore, exploiting the identities 



(CTI ■ r23)(o-2 • ^23) = -|^{5'23 + 0-23, CTls} 
„2 



(0-2 • ri3)(cri • ri3) = -^{o'23, -5*13 + CTis} , (43) 
6 

we can rewrite the Vd term in a form that has again 
the same spin-isospin structure as the anticommutator 
contribution of the UIX potential 

Vn{3 : 12) - ^{ri3, T23}[{ai3, a23}l^^^(ri3, r23)+ 



{'5'l3, Cr23}^J3^(''l3, ?'23) + 
Wl3, S23}Vo^ {ri3,r23)] , 



(44) 



VD^{ri3, r23) = y{ri3)zo{r23) + 2o(ri3)>'(?-23) 
^D^(^i3,?'23) = za{ri3)T{r23) 
VD^{ri3,r23) = T(ri3)zo{r23) ■ 



(45) 



In conclusion, including Vd amounts to properly adding 
the above radial functions to those already appearing in 
^3. 

The Ve term of TM' is completely equivalent to Vr (see 
Eq. p6p ). This allowed us to use the results of Ref. [2l| 
for the diagrams of Fig. 31 The same holds true for the 
chiral contact term Ve in PNM, as {Tij)pNM = 1, while 
in SNM the calculation of Ve requires the evaluation of 
the diagrams of Fig. [31 

The expression of diagram ([31 a) is 



(E1«)=Y^'EE / rf'r-i2dVi3^^,,i3^.'',' 



23" 



ex p 



^^,M2C(^B(3:12)Oy. (46) 

As pointed out in Ref. [2l| , integrating and tracing over 
the radial and spin-isospin variables of particle 3 leads to 
the appearance of an effective density dependent interac- 
tion 



E 



Vlyy',12iP)Ol2^pT. j ^''^13^^„,13^^V.23 >< 

ex 

C3fF£;(3:12)), (47) 



such that 



VV'A2 



(p). (48) 



ex p 



The subscripts xy label exchange patterns at the ends of 
the generalized correlation lines. In particular, dd cor- 
respond stands for direct-direct, de for direct-exchange, 
ee for exchange- exchange, and cc for incomplete circular 
exchange. The matrix Ap is defined through [l^l 



(49) 



implying AP = 1, 3, 3, 9, 6, 18 for p = 1, 6. 

It turns out that the only nonvanishing term of the 
density dependent potential is 



%yy'A2ip)=Wo''pJ2 I d'n3Z^y,13Zx'y',: 

ex 

Zo{ri3)zQ{r23) ■ 
The contribution of diagram ([Sib) is given by 



(50) 



ex p.p' 



[vE{i:l2)\{0P„0t,}] 



(51) 
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From the above expression, it clearly follows that only 
T-type generalized correlation lines contribute. Hence 



,23 



ex 

Zx"y",i2Zoiri3)zo{r23) ■ (52) 

Diagram ([31c) does not contribute to Ve, while for di- 
agram (ISld), with three generalized operatorial correla- 
tion, we find 



(P) 



.2^ „2 

' 2 



„^ „' nil 



ex p.p' .p' 



Ve(3:12) , 

x"y",23^\ gi 1*^121^13' 



0?3{0?2,Of3} + Of3{0?2,0?3})J. (53) 

The calculation of the spin-isospin traces yields 



(P) 



ce 



j dVi2dVi3(-2Z[2^[3Zj3 + 



^12 ^13^23 



^12^13 ^23 



ri ryiJT lyiJT ry(7T ■ -| Q^l^tt ryt rytT ryC7T 

^^^231^12^13^23 + J^°';231'^12^13 ^23 + 

1 Q(:tGt rytT rya ryt i ly t^tct rprtT ryaT rytT i 

^^S231^12'^13'^23 ^^^2?>\^\2^\Z ^23"'" 

(\^2^tG rya rytT rytT i Q/^iiCT ryCTT yt ryt 

^^231^12^^13^23 ^S231^12 '^13^23" 

r>(:tta ryOT ytT rytT 1 rt (^ttt ytT ytT ytT i 

"S231^12 '^13'^23 ^^S231'^12'^13'^23"^ 

-\QC'^tt ryt ytT rytT i 'XQt'^tt rytT ryt ryt i 

^-0^231^12^13^^23 J-^S231^12^13^23"r 
18^2^3*^ ^^2^13^23) ^0 (ri3 )^0 (r23 ) . 



(54) 



The matrices ^231 ^ depending on the angles formed by the 
vectors ri, Vi and r3, are defined in Ref. [i^. Following 
Ref. [l^l we have considered only the direct term of the 
generalized operatorial correlations. As a consequence, 



lyv 

dd,ij ' 



in the previous equation Zfj 

In order to find the optimal values of the variational pa- 
rameters, we have employed a procedure similar to sim- 
ulated annealing, the details of which are explained in 
Ref. 

The authors of Ref. 



26| constrained the difference be- 



tween the Pandharipande-Bethe (PB) and the Jackson- 
Feenberg ( JF) kinetic energies to be less than 10 % of the 
Fermi Energy Tp and required the sum rule involving the 
scalar two-body distribution function, g''(r 12), to be ful- 
filled with a precision of 3 %. In variational calculations 
of SNM we have imposed the further condition, firstly 
considered in Ref. |23j , that the sum rule of the isospin 
component of the two-body distribution function 



p / d^l25'^(''l2) 



-3, 



(55) 



be also satisfied to the same accuracy. 

Using also the sum rules for the spin and spin-isospin 
two body distribution functions leads to a sizable increase 



of the variational energies, which turn out to be much 
higher than those obtained releasing the additional con- 
straints, as well as the AFDMC results. The same pat- 
tern is observed in the results of variational calculations 
not including TNF. For this reason, we have enforced the 
fulfillment of the sum rules for g'^{ri2) and .g^(ri2) only. 

For potentials other than UIX, it turns out that the 
variational energies of PNM resulting from our optimiza- 
tion procedure are lower than the AFDMC values at 
p > Po . By carefully analyzing the contributions of the 
cluster expansion diagrams, we realized that the value of 
diagram ([3h) was unnaturally large. In particular, we 
have found that a small change in the variational pa- 
rameters leads to a huge variation of the value of the 
diagram. Moreover, the minimum of the energy in pa- 
rameter space was reached in a region where the kinetic 
energy difference was very close to the allowed limit. 

To cure this pathology, we have constrained the PB-JF 
kinetic energy difference to be less than 1 MeV, regard- 
less of density. The variational energies obtained impos- 
ing this new constraint are always larger than the corre- 
sponding AFDMC values and the value of diagram ([3lb) 
is brought under control. For the sake of consistency, 
the same constraint on the kinetic energies has been also 
applied to SNM. 



C. AFDMC calculations 

We have computed the EoS of PNM using the AFDMC 
approach 4l| with the TM' and NNLOL chiral potentials 
combined with the Argonne v'g NN interaction. 

An efficient procedure to perform AFDMC calculations 
with three-body potentials is described in Ref. [i^l ■ Since 
V3 is equivalent to the anticommutator term of the UIX 
model (while the commutator, V4, is zero in PNM), and 
in PNM the Ve terms of both the TM' and NNLOL po- 
tentials do not show any formal difference with respect to 
the repulsive term of UIX, the inclusion of these terms re- 
duces to a replacement of constants and radial functions. 
The authors of Ref. [42[ also described how to handle the 
Vi for the TM model, and no further difficulties arise in 
the case of the NNLOL potential. 

As the Vd term has never been encompassed in 
AFDMC, it is worthwhile showing how the calculation 
of this term reduces to a matrix multiplication. The ex- 
pectation value of Vd is given by 



(Vd) = J2 [^^(* • •?^) + ■■ «fc) + Voik : ij)] (56) 

i<j<k 

with Vd(* ■ jk) — Vd(* ■ kj) (otherwise all six permuta- 
tions need to be summed). Thanks to this property one 
can write 



{Vd)= VD{3:^k). 

i<k^j 



(57) 
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It is possible to write Vd(j : ik) of Eq. in terms of 
cartesian components operators 



(58) 



where 

















= T(r,,)(3ff/,^.-<5"^). 


(59) 



The anticommutation relation {crfiCr^} — 26"^ makes 
the expectation value of Vd a sum of 3iV x 3A^ matrix 
multiplications 

(Vd) = 2 ^ iXai-pjZpj-Sk + Zai-pjYfij-Sk + 
i<k,j 

Tai;PjZ[jj-Sk + Z ai-pjTpj-Sk)^? ^k 

- 2 ^({r, Z} + {T, Z})^,,sk <jT4 , (60) 

i<k 

analogous to those of Ref . [i^l ■ In order to compute the 
expectation value of Vd the former expression has been 
added to the cartesian matrices associated with the two- 
body potential. 

Following Ref. [H], we simulated PNM with A = 66 
neutrons in a periodic box, as described in Refs. [isl-lil]. 
using the fixed-phase approximation. For 66 neutrons 
finite-size effects on the kinetic energy have been found 
to be small, as its value is very close to the thermody- 
namic limit. Moreover, as shown in Ref. [43 |. the energy 
per particle obtained with 66 neutrons imposing the Pe- 
riodic Box Condition (PBC) differs by no more than 2% 
from the asymptotic value calculated with Twist Aver- 
aged Boundary Conditions (TABC). 

Finite-size effects are expected to be larger when the 
density is bigger, as the dimension of the box decreases. 
In order to check the validity of our calculations, at 
p = 0.48 fm~^ we have repeated the calculation with 
114 neutrons. For all potentials, no significant deviations 
from the energy per particle obtained with 66 neutrons 
have been observed. 
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IV. NUCLEAR MATTER EOS 
A. TM' potential 



Figure 5. Equation of state of PNM obtained using the 
AFDMC (triangles) and FHNC/SOC (solid lines) approaches 
with the TM' plus v'g hamiltonian. 



The results of Fig. [Sj showing the density dependence 
of the energy per nucleon in PNM, indicate that, once 
the new constraint on the difference between PB and 
JF kinetic energies is imposed, the agreement between 
FHNC/SOC (solid line) and AFDMC (triangles) resuffs 
is very good. 



The most striking feature of the results displayed in 
Fig. [S] is that, despite the parameters of the three body 
potentials being different, all SNM EoS obtained from 
the TM' potential turn out to be very close to each other. 
This is probably due to the fact that these potentials are 
designed to reproduce not only the binding energies of 
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and ^Hc, but also the n-d doublet scattering length 

aud- 
it is remarkable that although the parameters of TM' 
potentials were not adjusted to reproduce nuclear matter 
properties, the EoS saturates at densities only slightly 
lower than po = O-lSfm^"^, and the compressibilities are 
in agreement with the experimental value K sa 240 MeV. 
On the other hand, the binding energies are larger than 
the empirical value Eq — —16 MeV and rather close to 
the one obtained from the UIX potential, ~ 10 MeV [26l |. 
The numerical values of all these quantities are listed in 
Table El 



100 




0.1 0.2 0.3 0.4 0.5 

P (fm-^) 

Figure 6. Equation of state of SNM resulting from 
FHNC/SOC variational calculations with the TM' plus Ug 
hamiltonian. 



Table V. Saturation density, binding energy per particle and 
compressibility of SNM corresponding to the TM' EoS dis- 
played in Fig. |6l 





TM'i 


TMa 


TM3 


po (fni-3) 


0.12 


0.13 


0.14 


Eo (MeV) 


-9.0 


-8.8 


-9.4 


K (MeV) 


266 


243 


249 



B. NNLOL chiral potentials 

The results displayed in Fig. [7] show that, as in the 
case of the TM' potentials, the EoS of PNM computed 
within the AFDMC and FHNC/SOC schemes are very 
close to each other over the entire density range. 

The EoS of Fig. [7] are softer than those obtained from 
both the TM' (compare to Fig. [5]), and UIX (sse, e.g. 
Fig. 12 of Ref. [2al) potentials. This is due to the ambi- 
guity in the term Ve, discussed in Section UlI Al 



In the NNLOL2, NNLOL3, and NNLOL4 models the 
constant ce is negative. Therefore, the contribution of 
Ve is attractive, making the EoS very soft. When Ve is 
repulsive (i.e. ce is positive), as in the NNLOLi poten- 
tial, its contribution is very small and the resulting EoS, 
while being stifFer than those corresponding to the other 
NNLOL potentials, remains very soft. 

The recent astrophysical data of Ref. [lO] suggest 
that the EoS of PNM be at least as stiff as the one ob- 
tained with a readjusted version of the effective density- 
dependent potential of Lagaris and Pandharipande in 
combination with the Argonne v'q two-body interaction 
[45I I. Therefore, the EoS resulting from chiral NNLOL 
potentials are not likely to describe a neutron star of 
mass around 2Mq. 

The SNM EoS corresponding to the NNLOL potentials 
are displayed in Fig. [S] The fact that the NNLOL4 po- 
tential provides the stiffest EoS, while in PNM provided 
the softest, is not surprising. As discussed in Section 
nil Al when the contact term is attractive in PNM, it is 
repulsive in SNM, and viceversa. 

The results listed in Table IVII show that none of the 
chiral NNLOL potentials fulfills the empirical constraints 
on the SNM EoS. All potentials overestimate the satura- 
tion density, while the compressibility is compatible with 
the empirical value only for the NNLOL2 and NNLOL3 
models. As for the binding energies, they are closer to 
the experimental value than those obtained using both 
the UIX and TM' potentials. 

As a final remark, it has to be noticed that using the 
scalar repulsive term instead of provides more 
repulsion, resulting a stiffer EoS. As stressed in Section 
nil A| this issue needs to be addressed, taking into ac- 
count all terms that become equivalent in the limit of 
infinite cutoff only. Moreover, since the discrepancies 
among these terms are of the same order as the NNNLO 
term of the chiral expansion, other contact terms have to 
be included f46| . 

Table VI. Saturation density, the binding energy per particle, 
and the compressibility related to the NNLOL Eos displayed 
in Fig. [8] 





NNLOLi 


NNLOL2 


NNLOL3 NNLOL4 


po (fm-3) 


0.21 


0.20 


0.19 0.17 


Eo (MeV) 


-15.2 


-14.6 


-14.6 -12.9 


K (MeV) 


198 


252 


220 310 



V. CONCLUSIONS 

A new generation of chiral inspired three-nucleon po- 
tentials in coordinate space, suitable for carrying out nu- 
clear matter calculations, is now available. We have car- 
ried out a comparative analysis of the EoS of PNM and 
SNM obtained using the different parametrizations of the 
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Figure 7. Same as Fig. [5] but for NNLOL plus Vg hamiltonian. 
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Figure 8. Same as Fig. [6l but for NNLOL plus v'g hamiltonian. 



NNLOL potential, as well as the improved versions of the 
TM model discussed in Rcf. [11. 



The calculation of the SNM EoS has been been per- 
formed within the variational FHNC/SOC approach. In 
the case of PNM we have also used the AFDMC com- 
putational scheme, the results of which turn out to be 
in close agreement with the variational FHNC/SOC es- 
timates. 

Our analysis shows that the transformation from mo- 
mentum to coordinate space brings about a cutoff depen- 
dence, leading to sizable effects in nuclear matter. As 
discussed in Section IIII Al the contribution of the con- 
tact term, which in PNM would vanish in the A — ?> co 
limit, can not be fully determined fitting the low energy 
observablcs. Moreover, the NNN contact terms of the 
NNLOL2 and NNLOL3 models turn out to be attractive 
in PNM, leading to a strong softening of the EoS. 

An illustrative example of the uncertainty associated 
with the local form of the NNN contact term is provided 
by the results of Fig. [8] and Table ED The NNLOL4 
model largely overestimates the empirical value of the 
compressibility modulus of SNM, thus yielding a stiff 
EoS. On the other hand, as pointed out in Section IIVBI 
it predicts a soft EoS of PNM. The impact of this is am- 
biguity is large, since compressibility is a most important 
property of the EoS. The recent discovery of a '--^ 2 M0 
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neutron star appears in fact to rule out dynamical models 
yielding a soft EoS of /3-stable matter. 

None of the considered three-nucleon potential models 
simultaneously explains the empirical equilibrium density 
and binding energy of SNM. However, among the differ- 
ent parametrization that we have analyzed, the NNLOL4 
and TM3 provide reasonable values of po- It has to be 
emphasized that this is a remarkable result, as, unlike the 
UIX model, these potential do not involve any parameter 
adjusted to reproduce po- 

In order to resolve the inconsistencies involved in the 
contact term, one should include all contributions to this 
term arising from the chiral expansion at NNLO. More- 
over, as pointed out by the authors of Rcf. [i^ , due to the 
choice of the regulator function (see Eq.®), a fully con- 
sistent treatment should also take into account NNNLO 
contact contributions. 

As a final remark, it must be mentioned that the TM' 
and NNLOL potentials discussed in this paper can be 



used to obtain two-body density-dependent effective in- 
teraction within the formalism developed in Ref . [IBl • At 
present; this is the only approach allowing for the inclu- 
sion of three-nucleon potentials involving a term of the 
form of Vi of Eq.© in AFDMC calculations of SNM. 
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